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Data collected by the DO detector at a pp center-of-mass energy of 1.96 TeV at the Fermilab 
Tevatron Collider have been used to search for pair production of the lightest supersymmetric 
partner of the top quark decaying into blv. The search is performed in the II' = e/i and pLpL final 
states. No evidence for this process has been found in data samples of approximately 400 pb~^. 
The domain in the [M {ti) , M {i>)] plane excluded at the 95% C.L. is substantially extended by this 
search. 

PACS numbers: 14.80.Ly; 12.60.Jv 



Supersymmetric theories [l| predict the existence of 
a scalar partner for each standard model fermion. Be- 
cause of the large mass of the standard model top quark, 
the mixing between its chiral supersymmetric partners is 
the largest among all squarks; therefore the lightest su- 
persymmetric partner of the top quark, ti (stop), might 
be the lightest squark. If the ti hlv decay channel is 
kinematically accessible, it will be dominant y as long as 
the ii — > bxi and ii tx\ channels are kinematically 
closed, where Xi ^'Hd Xi are the lightest chargino and 
neutralino, respectively. In this letter we present a search 
for stop pair production in pp collisions at 1.96 TeV with 
the DO detector, where a virtual chargino decays into 
a lepton and a sneutrino, and where the sneutrino v, 
considered to be the next lightest supersymmetric par- 
ticle, decays into a neutrino and the lightest neutralino 
x\] in pp collisions, stop pairs are dominantly produced 
via the strong interaction in quark antiquark annihila- 
tion and gluon fusion. We use the Minimal Supersym- 
metric Standard Model (MSSM) as the phenomenologi- 
cal framework for this search. We assume the branching 
ratio Br{xi ~* — 1 with equal sharing among all 
lepton flavors, and we consider only cases where (. = e, 
For stop pair production, we consider hh iCiyVxiXi fi- 
nal states with ££' = e^/i^ and ££' = {en and 
/i/i channels) ; the signal topology consists of two isolated 
leptons, missing transverse energy (-^t), and jets. DO has 
also searched for scalar top in the charm jet final state 

i- 

The DO detector Q comprises a central tracking sys- 
tem surrounded by a liquid-argon sampling calorimeter 
and a system of muon detectors. Charged particles are re- 
constructed using a multi-layer silicon detector and eight 
double layers of scintillating fibers in a 2 T magnetic field 
produced by a superconducting solenoid. The calorime- 
ter provides hermetic coverage up to pseudo-rapidities 
\ri\ ~ 4 (where 7y=-log(tan(6'/2), and where 6 is the po- 
lar angle with respect to the proton beam direction) in 
a semi-pro jective tower geometry with longitudinal seg- 
mentation. After passing through the calorimeter, muons 
are detected in the muon detector comprising three lay- 
ers of tracking detectors and scintillation counters located 
inside and outside of 1.8 T iron toroids. Events contain- 
ing electrons or muons are selected for off-line analysis 
by a trigger system. A set of dilepton triggers is used to 
tag the presence of electrons and muons based on their 
energy deposit in the calorimeter, hits in the muon de- 



tectors, and tracks in the tracking system. 

Three-body decays of the ti are simulated us- 
ing COMPHEP 5] and PYTHIA [6] for generation and 
hadronization respectively. Standard model background 
processes are simulated using the pythia and ALPGEN 
,7;] Monte Carlo (MC) generators. These MC samples are 
generated using the CTEQSl [sI parton distribution func- 
tions (PDF); they are normalized using next-to- leading 
order cross sections 9]. All generated events are passed 
through the full simulation of the detector geometry and 
response based on GEANT [10] . MC events are then recon- 
structed and analyzed with the same programs as used 
for the data. 

Muons are reconstructed by finding tracks pointing 
to hit patterns in the muon system. Non-isolated 
muons are rejected by requiring the sum of the trans- 
verse momenta (pt) of tracks inside a cone with ATZ — 
^ (A(/))2 -I- (Ayy)^ ~ 0.5 (where <^ is the azimuthal angle) 
around the muon direction, and the calorimeter energy 
in an annulus of size 0.1 < AT?, < 0.4 around the muon to 
be less than 4 GeV/c and 4 GeV. Isolated electrons are 
selected based on their characteristic energy deposition 
in the calorimeter, their fraction of deposited energy in 
the electromagnetic portion of the calorimeter and their 
transverse shower profile inside a cone of radius AT?. = 0.4 
around the direction of the electron; furthermore, it is re- 
quired that a track points to the energy deposition in the 
calorimeter and that its momentum and the calorimeter 
energy are consistent with the same electron energy; an 
"electron-likelihood" is defined as a variable combining 
information from the energy deposition in the calorime- 
ter and the associated track. Backgrounds from jets and 
photon conversions are further suppressed by requiring 
the tracks associated with the muons and electrons to 
each have at least one hit in the silicon detector. Jets are 
reconstructed from the energy deposition in calorimeter 
towers using the Run II cone algorithm ll| with radius 
AT?=0.5, and corrected for the jet energy scale (JES) 
[l^ : in this search, jets are considered with pT > 15 
GeV/c. The I^t is defined as the energy imbalance of all 
calorimeter cells in the plane transverse to the beam di- 
rection, and is corrected for the JES, the electromagnetic 
energy scale, and reconstructed muons. All efficiencies 
are measured with data \2M- 

In both e/i and fifj, channels, the signal points 
[M{ii),M{i>)] = (110,80) GeV/c2 and (145,50) GeV/c^, 
respectively referred as "soft" (point A) and "hard" 
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(point B) signals, have been used to optimize the selec- 
tion of signals of different kinematics because of different 
Am = M{ii) — M{v). The choice of these points was 
also motivated by the sensitivity of the DO search during 
Run I [l^ . The main background processes imitating the 
signal topology are Z/^*, WW, production, and mul- 
tijet background. All but the latter are estimated with 
MC simulation. The multijet background is estimated 
from data. In the e/i channel, two samples each domi- 
nated by a different multijet background are obtained by 
inverting the muon isolation requirements, and by invert- 
ing the cut on the electron-likelihood; in the /i/i channel, 
such a sample is obtained by selecting same-sign muon 
events. Factors normalizing each sample to the selection 
sample are also obtained from data, and applied to the 
background samples to obtain the multijet background 
estimation, this, at an early stage of the selection. 

For the e/i channel, the integrated luminosity [31 of 
the data sample is (428 ± 28) pb""'^. The preselection 
is concluded by requiring the transverse momenta of the 
electron and muon (see Fig. [TJa) and (b)) to be greater 
than 10 and 8 GeV/c, respectively. In this final state, the 
data are dominated by the multijet and Z /^* — > rr back- 
grounds. In these processes, poorly reconstructed leptons 
are correlated with IpT, giving rise to higher event pop- 
ulations at high and low values of the azimuthal angular 
difference between the leptons and the ^t, a low value 
of the angular difference for one lepton being correlated 
with a high value of the other. Taking advantage of a 
higher background contribution at low values of angular 
distributions, we require 

A(/)(/i,^T) > 0.4 , A0(e,^T) > 0.4. (Emu 1) 

We require If^r to be greater than 15 GeV to reduce 
contribution of both the multijet and Z/^* tt back- 
grounds. To reject multijet events in which leptons are 
associated with a jet, we require the two leptons to be at 
a ATZ distance greater than 0.5 from any reconstructed 
jet. To further reduce the multijet contribution, we re- 
quire the z component of the origin of the highest trans- 
verse momentum muon track to be within four standard 
deviations a from the z component of the primary vertex: 

> 15 GeV 
A7^[(e,/i),jet] > 0.5 
|z(/i) - z(p.v.)| < 4cr. (Emu 2) 

To reduce the Z/j* — > tt background, we cut on 
low values of the transverse mass of the muon and I^t 
{MxifJ,, ]^t), see Fig. [TJc)). To further reduce this back- 
ground, we make use of the correlation between the an- 
gular differences A(j){iJ,, ]^t) and A0(e, .^t), and require 
their sum (see Fig. [Ijd)) to be greater than 2.9: 

MT(/i,^T) > 15 GeV/c2 
/\(j){^i,$T) + A(t>le,:$T) > 2.9. (Emu 3) 



The contributions of different backgrounds, and the ex- 
pected numbers of signal and observed data events in the 

final state at different selection levels are summarized 
in Table [H After all selections, the WW (dominating the 
diboson contribution) and tt contributions are the domi- 
nant backgrounds. To separate soft signals such as point 
A from these backgrounds, we consider the variable St 
defined as the scalar sum of the transverse momentum of 
the muon, the electron, and the (see Fig. [He)). To 
separate hard signals such as point B from background 
contributions, we consider the variable Ht defined as the 
scalar sum of the transverse momentum of all jets (see 
Fig. [Hf)). Rather than cutting on these two variables, 
the Ht and St spectra predicted for signal and back- 
ground are compared with the observed spectra in twelve 
[St, Ht] bins (see Table HI)) when extracting limits on the 
signal cross section, thus allowing a separation of signals 
of different kinematics from the WW and tt backgrounds. 

For the fifi channel, the integrated luminosity [l5l | of 
the data sample is (395 ± 26) pb~^. The selection of the 
signal in this final state is more challenging because of 
the strongly dominating Z/^* background. The 

preselection is concluded by requiring the transverse mo- 
menta of the two highest transverse momenta opposite- 
sign muons to be greater than 8 and 6 GeV/c. While the 
signal is characterized by the presence of jets originat- 
ing from the hadronization of b quarks, the Z/j* — > fifi 
background owes the presence of jets to initial state radi- 
ation gluons which hadronize into softer jets, resulting in 
a lower multiplicity of jets; the latter is also valid for soft 
signals such as point A. To keep sensitivity to soft sig- 
nals while rejecting substantial background, we require 
at least one jet: 

7V(jets) > 1. (Dimu 1) 

To further remove Z/j* — )■ /i/i background events, where 
poorly reconstructed muons correlate with the J^t, we 
require the ]^t to be greater than the contour shown 
on Fig. m^a), using a cut parametrized by the following 
equation: 

/ GeV > 20+ |A0(/ii,^T) - 1.55|9-2, (Dimu 2) 

where fii is the highest transverse momentum muon. To 
augment the search sensitivity in this channel, we take 
advantage of the presence of jets originating from the 
fragmentation of long-lived b quarks in the signal. An 
algorithm based on the lifetime of hadrons calculates the 
probability V for the tracks of a jet to originate from 
the primary interaction point [l^ . This b jet tagging 
probability is constructed such that its distribution is 
uniform for light-flavor jets while peaking at zero for 
heavy-flavor jets which have a vertex significantly dis- 
placed from the primary vertex (Fig. [l^b)). Considering 
the highest transverse energy jet, we require 

P(jet) < 1%. (Dimu 3) 
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A cut on the dimuon invariant mass (Fig. [H^c)) in the 
vicinity of the Z boson resonance only at low Jf^T (Fig- 
[21^d)) further suppresses the Zj^* /x/i background 
while preserving the signal: 

Af (/i, /i) ^ [75, 120] GeV/c2 for < 50 GeV. (Dimu 4) 



Table IIIII summarizes the different stages of the signal 
selection in the /i/i channel. The tt background domi- 
nates after the selection cuts; five Ht bins are considered 
(see Table WV} to separate various signal points from this 
background. 

The expected numbers of background and signal events 
depend on several measurements and parametrizations 
which each introduce a systematic uncertainty: lep- 
ton identification and reconstruction efficiency [(2.6-7)%] 
13 1, trigger efficiency [(3.5-5)%] [l^, luminosity [6.1%] 
15 1, multijet background modeling [10%], JES [(4-22)%] 
12l ]. jet identification and reconstruction efficiency and 
resolution [(4-16)%] [H, b jet tagging [(1-11)%] [11, 



PDF uncertainty affecting the signal efficiency [10%] [17|. 

After applying all selection cuts for e/i and fi/j, data 
sets, no evidence for ti production is observed. We 
combine the number of expected signal and background 
events and their corresponding uncertainty, and the num- 
ber of observed events in data from the twelve bins of the 
en selection ( Table lll| and the five bins of the fifi selection 
(Table lIVp to calculate upper-limit cross sections for sig- 
nal production at the 95% C.L. for various signal points 
using the modified frequentist approach [3]. In this cal- 
culation, correlated uncertainties are taken into account; 
no overlap is expected nor observed between the two sam- 
ples. Regions for which the calculated cross section upper 
limit is smaller than the theoretical one are excluded at 
95% C.L. Figure [3] shows the excluded region as a func- 
tion of the scalar top quark and sneutrino masses, for 
nominal (solid line) and for both minimal and maximal 
(band surrounding the line) values of the titi produc- 
tion cross section; the latter variation corresponds to the 
PDF uncertainty for the signal cross section, quadrati- 
cally added to the 2/^^ and /i, /2 renormalization scale 
variations of the titi cross section. Although the num- 
bers of expected and observed events are similar (Tables 
|T] and , their distribution across the bins (Tables |TT] 
and IIV|1 causes the expected cross section limit to be 
lower than the observed one. For minimal values of the 
production cross section, the search in the e/i final state 
individually excludes a stop mass of 176 GcV/c^ for a 
sneutrino mass of 60 GeV/c^, and a sneutrino mass of 
97 GeV/c^ for a stop mass of 130 GeV/c^; the search in 
the fifj, final state, once combined with the efj, final state, 
extends the final sensitivity by approximately 10 GeV/c^ 
for small and large mass differences. 

In summary, we have searched for the lightest scalar 
top quark decaying into Miy] events with an electron and 
a muon, and two muons have been considered for this 



search. No evidence for the lightest stop is observed 
in these decays, leading to a 95% C.L. exclusion in the 
[M{ii), M{i>)] plane. The largest stop mass excluded is 
186 GeV/c^ for a sneutrino mass of 71 GeV/c^, and the 
largest sneutrino mass excluded is 107 GeV/c^ for a stop 
mass of 145 GeV/c^; these mass limits are obtained with 
the most conservative theoretical production cross sec- 
tion, taking into account the PDF uncertainty and the 
variation of the renormalization scale. This is the most 
sensitive search for stop decaying into Mi> to date. 
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TABLE I: en channel. Expected numbers of events in various background and signal channels, and number of observed events 
in data, at various selection levels. Statistical as well as systematic uncertainties from the JES correction are shown for the 
total background and signal. 



Scloctiini 


Multijot 


Background 

^/v - 


contributions 

a if 


Dibosou 


Total 
Ba.ckgrinuKl 


Data 


s: 

Point .1 


ignal 

Point B 


Preselection 
Emu 1 
Emu 2 
Emu 3 


304.5 
194.4 
8.6 
5.9 


286.7 
115.4 
20.0 
3.6 


12.4 
10.4 
9.1 
7.4 


28.6 
25.3 
21.2 
20.2 


632.3 ± 19.5115;" 

345.4 ± 15.0l!!;7 
58.9±3.8l2:2 
37.1 ± 2.7l°| 


596 
329 
52 
34 


65.9 ± 2.4i;;;'^ 
54.1 ±2.21!',;;; 
31.6± 1.7lo:o 
26.0 ± l.f>tol 


26.6 ± 0.7l!!;^ 
22.7±0.7in 
19.0±0.6lo:? 
17.3 ± 0.6±°| 



TABLE IL ejj, channel. Expected numbers of events for total background, signal points A and B, and number of observed 
events in data, in the twelve [St, Ht\ bins. Statistical and JES uncertainties are added in quadrature for the total background 
and signal points. 



Total Signal 



Bin 


background 


Data 


Point A 


Point B 


St e [0, 70[ GeV, = 


2.6 ± 1.1 


1 


7.3 ±1.0 


0.0 ±0.0 


St G [70, 120[ GeV, Ht=Q 


9.2 ± 1.2 


14 


4.8 ± 0.7 


0.2 ±0.1 


St € [120, ...[ GeV, i?T = 


7.7 ±0.7 


5 


0.8 ±0.3 


1.8 ±0.2 


St € [0, 70[ GeV, Ht £ ]0, 60] 


1.9 ±0.7 


2 


5.2 ±0.7 


0.0 ±0.0 


St € [70, 120[ GeV, Ht £ ]0, 60] 


3.6 ±1.2 


4 


5.3 ±0.8 


1.2 ±0.2 


St e [120, ...[ GeV, Ht G ]0,60] 


3.0 ± 0.4 


2 


0.6 ±0.3 


6.3 ±0.5 


St G [0, 70[ GcV, Ht G ]60, 120] 


0.4 ±0.6 





0.6 ±0.3 


0.0 ±0.0 


St G [70, 120[ GoV, Ht G ]60, 120] 


0.7 ±0.2 


1 


1.2 ±0.3 


1.3 ±0.2 


St G [120, ...[ GeV, Ht G ]60, 120] 


3.6 ±0.8 


2 


0.1 ±0.1 


4.3 ±0.3 


St G [0, 70[ GeV, Ht G ]120, ...[ 


0.0 ±0.0 





0.0 ±0.0 


0.0 ±0.0 


St G [70, 120[ GoV, Ht G ]120, ...[ 


0.8 ±0.6 


1 


0.0 ±0.0 


0.4 ±0.1 


St G [120,...[ GeV, Ht G ]120,...[ 


3.7 ± 1.1 


2 


0.1 ±0.1 


1.7 ±0.3 



TABLE III: fifi channel. Expected rmrnbors of events in various background and signal channels, and rmnibcr of observed 
events ill (lata, at various soloctiori levels. Statistical as well as systeiiiatic: uncertainties from the JES corre(;tioii are shown for 
the total hackgrouud and signal. 



Selection 


Multijet 


Background contributions 
T(l,25) Z/-f*^U ti 


WW 


Total 
Background 


Data 


Signal 

Point A Point B 


Preselection 
Dimu 1 
Dimu 2 
Dimu 3 
Dimu 4 


3607.6 
682.1 
41.8 
0.0 
0.0 


973.1 
80.8 
0.4 
0.0 
0.0 


23781.7 
3894.9 
155.7 
6.1 
0.1 


5.1 
5.1 
4.7 
2.6 
2.3 


9.6 
1.5 
1.1 
0.0 
0.0 


28377.1 ±348lo:o 
4664.4 ± 97l^g 
203.7 ± Stll 
8.7± 1.6tJ? 
2.9 ± OAtoA 


28733 
4337 
213 
4 
1 


9.8 ± OAto^^ 
8.8 ± 0.4t°;l 
7.5 ±0.31°:? 
3.5 ± 0.2l° g 
3.1 ± 0.2l^ o 


4Ll±1.5l»:g 
24.2 ± 1.111:1 
12.9 ±0.811-3 
3.4 ±0.41°:^ 
3.3±0.4l°:* 



TABLE IV: nfj, channel. Expected numbers of events for total background, signal points A and S, and number of observed 
events in data, in the 5 Ht bins. Statistical and JES uncertainties are added in quadrature for the total background and signal 
points. 



Bin 


Total 
background 


Data 


Point A 


Signal 


Point B 


Ht G jo, 40] GeV 


0.11 ±0.0 





2.0 ±0.3 




0.5 ±0.1 


Ht G ]40, 80] GeV 


0.89 ± 0.4 





1.1 ±0.3 




1.0 ±0.1 


Ht G ]80, 120] GeV 


0.75 ± 0.0 





0.2 ±0.1 




0.8 ±0.1 


Ht G ]120, 160] GeV 


0.56 ±0.0 


1 


0.0 ±0.0 




0.4 ±0.1 


Ht G ]160, ...[ GeV 


0.57 ± 0.0 





0.0 ±0.0 




0.4 ±0.1 
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FIG. 1: e/i channel. Distributions of the transverse momenta of the electron (a) and of the muon (b) after preselection cuts; 
(c) the transverse mass Mt{^j.,IPt) after preselection cuts and IpT > 15 GeV and Ai?[(e,^),jet] > 0.5; (d) the angular sum 
A(I>{^,1^t) + A(l>{e,l^T) after the cut (Emu 2); (e) St and (f) Ht distributions after the cut (Emu 3). 
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FIG. 2: /i/i channel, (a) A</)(/xi, _^t) versus in simulated Z/j* fijj, events; the contour of the cut (Dimu 2) is shown by 
the solid line. Distributions of the 6 jet tagging probability ^(jet) (b), the invariant mass of the two most energetic muons (c), 
and (d) after preselection cuts. 
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DO Run II observed limit: 

(e,|.L) 428 pb'^ and (il^l) 395 pb'^ 




60 80 100120140160180 200 220 

Scalar top quark mass [GeV/c^] 



FIG. 3: For the nominal production cross section, the 95% C.L. excluded regions in the [M (ti) , M {{>)] plane for the observed 
(fuU curve) and the average expected (dashed curve) limits are shown; the band surrounding the observed limit represents the 
lower and upper bounds of the signal cross-section variation. The regions excluded by DO during Run I [ij and by LEP 
are also shown. 
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